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ABSTRACT 

The discovery (Liang & Kargatis 1996), that the peak energy of time-resolved spectra of gamma-ray burst 
(GRB) pulses decays exponentially with fluence, is analytically shown to imply that the time-integrated photon 
number spectrum of a pulse should have a unique shape, given by an underlying E~ x behavior. We also show 
that the asymptotic low energy normalization of the time-integrated spectrum is equal to the exponential decay 
constant. We study analytically how this general behavior is modified in more realistic situations and show 
that diversity is then introduced in the properties of time-integrated GRB pulse spectra. We argue that further 
diversity will occur in time-integrated multi-pulse (complex) GRB spectra. The total energy received per cm 2 is 
approximately the decay constant times the maximum peak energy of the pulse. Our analytical results connect 
the properties of the time-integrated pulse spectrum with those of the time-resolved spectra, and can thus be used 
when studying observed GRB pulse spectra. We illustrate with the bright burst GRB 910807 and comment on 
GRB 910525 and GRB 921207. 

Subject headings: gamma rays: bursts 


1. INTRODUCTION 

The time-integrated and time-resolved spectra of gamma-ray 
bursts (GRB) differ. The general appearance of time-integrated 
GRB spectra is that of a smoothly broken power law. Time- 
integrated photon number spectra from missions prior to the 
Compton Gamma-Ray Observatory (CGRO) with its Burst and 
Transient Source Experiment (BATSE) were often fitted with 
an “optically thin thermal bremsstrahlung” (OTTB) spectrum, 
Ne{E) oc E~ 1 e~ E ^ E °, combined with a power law at larger 
energies. Band et al. (1993) showed that nearly all time- 
integrated BATSE spectra were well fitted by the “Band” func¬ 
tion consisting of two power laws, JVe(E) oc E a at low ener¬ 
gies and Ne(E) oc at high energies, smoothly joined by 
an exponential roll-over characterized by a constant Eq. The 
Band function index, a, shows a broad peak between -1.25 and 
-0.25 (Band et al. 1993), instead of the earlier claimed univer¬ 
sal value of -1, /3 clusters fairly narrowly around -2.1 (Preece 
et al. 1997), while Eq has values between 100 keV and 1 MeV, 
but mainly clustering around 100-200 keV (Band et al. 1993). 

Time-resolved spectra, even within single pulses, show 
strong time evolution, and are, in general, harder than time- 
integrated spectra (e.g., Pendleton et al. 1994a, Ford etal. 1995, 
Liang & Kargatis 1996, hereafter LK96, Crider et al. 1997a, 
1997b). Most importantly, LK96 found that the peak energy, 
E’pk, of the time-resolved E 2 Me{E) spectra of a single pulse 
in a number of bright and smooth GRBs, decays exponentially 
as a function of photon fluence, i. e ., 

£pk(t) =^pk,maxe- 4 ' (t)/$0 . (1) 

where E/ p k,max is the maximum value of E p k within a pulse, 
$(f) is the photon fluence integrated from the time of -E p k,max, 
and (ho is the exponential decay constant. Furthermore, Crider 
et al. (1997a) found the instantaneous Band function index, cq, 
in general to be hard, > —1, reaching values as large as +1.6. 

The aim of this paper is to show the consequences of the 
time evolution of E p ^(t), for the time-integrated GRB spec¬ 
trum of a single pulse, as well as for a number of pulses, and 


to determine the meaning of >!>o. The following notation is 
used: N(t) (photons cm -2 s -1 ) is the photon number flux of 
the burst at time t, and Ne{E,1) (photons cm -2 s -1 keV -1 ) 
is the instantaneous spectrum of the burst at photon energy E 
(keV) at time t. Here, Ne(E, t)dE is the photon number flux 
at time t within an energy interval dE around E. Furthermore, 
A/e {E) = Ne(E, t) dt (photons cm -2 keV -1 ) is the time- 
integrated spectrum of the burst and E p k is the peak photon en¬ 
ergy of E 2 Ne(E, t) spectra. The total photon number fluence 
is $(f) = j p N{t')dt' (cm -2 ). Note that d$(t)/dt = N(t). 

In § 2, we show analytically that the exponential E p k- 
decay implies that the time-integrated spectrum of a single 
pulse should have a typical asymptotic low energy behavior of 
A/e (E) es $ 0 E~ l , i.e ., the same low energy power law slope 
as the OTTB spectrum (the “universal” a = — 1), and with the 
asymptotic low energy normalization given by <f>o- Modifica¬ 
tions to this behavior below E’pk.max are shown to be approx¬ 
imately given by an incomplete gamma function. In § 3, we 
discuss the implications for the overall spectra of GRBs and in 
§ 4, we illustrate our analytical findings by analyzing spectral 
data from the bright burst, GRB 910807. In particular, we show 
that the decay constant, $o, °f the spectral evolution is consis¬ 
tent with the asymptotic normalization of the time-integrated 
spectrum. Finally, in § 5, we discuss the effects of a time- 
varying spectral shape and illustrate using GRB 950525 and 
GRB 921207 


2. THE TIME-INTEGRATED SPECTRUM OF A SINGLE PULSE 

We present a scenario in which the shape of the time- 
integrated spectrum is a direct consequence, and therefore a 
probe of, the specific time-evolution of the instantaneous spec¬ 
tra. We assume the latter to be peaked in the ENe spectrum, 
i.e., have a peak energy E p around which the number of emit¬ 
ted photons dominate, which is equivalent to that the power law 
indices a > — 1 and (3 < —1. This is indeed the case in most 
studied pulses (see, e.g., Crider et al. 1997a), except at the very 
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end of the pulse decay, when slopes with a < —1 sometimes 
are observed. This does not necessarily contradict the follow¬ 
ing discussion, as it can be ascribed to the possible presence of 
a soft component (Preece et al. 1996) and the contributions of 
previous (possibly untriggered or unresolved) pulses. It should 
also be noted that the determination of a becomes more diffi¬ 
cult the closer the peak energy is to the low energy boundary of 
the observed spectral range. In any case, such a soft contribu¬ 
tion should be negligible compared to the total fluence for the 
discussion below to be useful. 

Integrating the instantaneous spectra over some time interval 
(t i to A) of the pulse, the time-integrated spectrum becomes 


a r E (E) 


J 2 N E (E,t)dt = 


-Epk ,mi n 

Ne{E, t(E pk ))—— dE pk , (2) 

W-C/pk 


' ^pk, 


where E p k,max — E p k (t — tl ) and Epk,min — E p k(f — A)* If 

Equation (1) is valid during the time interval being considered, 
we have 


dE pk (t) 

dt 


—E pk (t) 


m 

$0 


-E pk (t) 


N(E pk (t)) 

$0 


(3) 


where we introduce a new function, N(E pk (t)) = N(t). We 
note that the energy dependence of N E (E,t(E pk )) normally 
is on the ratio E/E pk , but that its spectral shape could be 
time dependent, adding additional dependence on t(E pk ) or 
E pk . We therefore define a function N e (E/E pk , E pk ), where 
N E (E,t) = N E (E/E pk ,E pk ) = N(E pk )f E (E/E pk ,E pk ), 
and where the function f E (E/E pk , E pk ) describing the 
instantaneous spectral shape is normalized to unity, 
f f E (E/E pk , E pk ) dE = 1. Then the time-integrated spec¬ 
trum, Equation (2), becomes 


fBpk.max ()P , 

Ms(E) = $ 0 / h(E/E pk , E pk )—?E. 

•/-Epk.min 


(4) 


Introducing f x (x,E pk ) = f E (E/E pk , E pk )dE/dx = 
E pk f E (E / E pk , E pk ), where x = E/E pk , and f x (x,E pk ) is 
normalized to unity, f f x (x, E pk ) dx = 1, we find that 


A4(E) = • 1(E), 

where the correction factor 


i(e) = r 

J Xrt 


*(E) 


n(E) 


f x (x,E/x) dx, 


(5) 

(6) 


is its asymptotic low energy normalization and the “universal” 
time-integrated a t i = — 1 is obtained. The only assumption 
made for the derivation is the empirical relation (1) shown so 
far to be valid for numerous burst pulses (LK96, Crider et al. 
1997a,c), but whose exact validity must be explored. 

Equation (1) means that a linear increase in fluence by equal 
steps of $0 photons cm -2 corresponds to a decrease of In E p in 
equal logarithmic steps. The instantaneous spectra are, roughly, 
dominated by the about $0 photons cm~ 2 in a logarithmic in¬ 
terval around E = E pk , and as a result EM e (E) is a constant 
= <1» 0 . As EM e (E) also is the time-integrated specific flux 
spectrum, it is clear that the photon spectrum, AT e (E) of a sin¬ 
gle pulse has a slope of —1. The OTTB-like spectrum is thus a 
result of the exponential decay of the peak energy versus photon 
fluence. 

Let us now study the behavior of 1(E) when we can neglect 
the time-dependent changes of the instantaneous spectral shape, 
i.e, when the E/x dependence of f x (x, E/x) in Equation (6) 
can be neglected. It has been shown for several bursts that the 
time-resolved spectra are well fitted by the Band function and 
therefore it is of interest to use the Band function in (6). How¬ 
ever, before discussing these results, we study two conceptu¬ 
ally simpler approximations, which are valuable in illustrating 
the resulting behavior of the time-integrated spectrum. Firstly, 
in the case of integrating over all energies, i.e., ,T ni j tl = 0 and 
x m ax = 00 , 1(E) = 1 (due to the normalization of f x (x)) 
and AT e (E) = $0 /E is an exact solution of (5). The shape 
of the instantaneous spectra does not affect the normalization 
and shape of the time-integrated spectrum as long as the instan¬ 
taneous EN E (E,t) spectra are peaked with the Band indices 
Oi, > — 1 and A < —1. Then the dominant contribution to the 
integral comes from x ~ 1 and I(E) is convergent. Secondly, 
if we consider finite integration limits in (6) and replace the in¬ 
stantaneous spectra with a Dirac delta function at E = E pk , 
i.e., x = 1 and f x (x) = 6(x — 1), we arrive at a similar result: 
1(E) = 1 for £pk,min < E < -Epk,max and 0 otherwise. 

Instantaneous spectra broader than a delta function will intro¬ 
duce deviations from the E -1 power law found in the previous 
case. We now consider the Band function in (6). Changing vari¬ 
ables to y = E/Eq and thus y = x(2 + a ;) the Band function 
becomes 

f [Aj-^e-y if y < (a-, - A); 

fy( y ) = < 

[ [A'\ 1 y f3i if y > («i - A). 

For f y to be a continuous function, we have 

jy = ^(“i - A)“ i_A e _( “ i_A) = jC( ai - A), (7) 

and the normalization, f f y dy = 1, gives 


and where x lnin (E) = E/E pk ,max and x max (E) = 
E /Epk,min- We show below that the E-dependence of 1(E) 
is rather weak between E p k, m in and E p k, m ax for the cases we 
consider. The first factor in (5), describing the over-all spec¬ 
tral shape, shows that the time-integrated spectrum of a single 
pulse, which has an exponentially decaying peak energy, is a 
power law with a time-integrated spectral index of at\ ~ — 1 
between E p k, m in and E p k, m ax normalized by the decay con¬ 
stant $ 0 - This result points out how the instantaneous and the 
time-integrated spectra are related. It shows that the exponen¬ 
tial decay behavior of the instantaneous spectra leaves a signa¬ 
ture on the time-integrated spectrum: the decay constant, Tq, 


A 


r(cq + l)[P(aj + 1, cq 
e -(ai-/3i) 

T(cq + 1)(A + 1 ) 



( 8 ) 


where P(a ; + 1 ,y) is the incomplete gamma function (see, 
e.g.. Press et al. 1992). We have three cases: 

(i) 2/max > Oi - A) and y min < (a ; - A) gives 


1(E) = r (°^ 4 + ^ [POi + 1, «i - A) ~ P(cn + 1, 2/min) 

+ r ( „?+i)(ali) ( ^^ ( “’-^ )fl+ 4 (9) 
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Fig. 1.— The correction factor, 1(E), given by Eq. (6) plotted as a function of E /-Epk.max for different values of the instantaneous 
o:j. The curves are for f x (x) being a Band function, but the principle behavior below E p k. m ax is fairly well described by Equation 
(10) alone, i.e, for f x (x) being an exponentially cut-off power law. The solid curves describe a situation with no low energy cut-off, 
-Epk.min = 0. Then Eq. (10) becomes 1(E) = [T(ai + 1)/^4][1 — P(a i + l,t/ m in)]- The dotted curves include such a cut-off, 
assuming £ p k,max/-Epk,min = 20. This corresponds to using the full Eq. (10). 


(ii )y max < («i - Pi) gives 


integrated spectrum approaches 


1(E) 


T(ai+1) 

A 


\P (oii ~\~ 1, 2/max) P{pt\ H - 2/min)] • 

( 10 ) 


This is the spectrum below (cq — 8 \) Eq , and it is independent 
of Pi. 

(iii) 2/min > (cq - Pi) gives 


1(E) 


rOi + i) 

A 


C'(aj Pi) , a+i 

r(ai + i)(A + i)^ max 



At 2 / max = (cq — Pi) cases (i) and (ii) are identical, and at 
2/min = (ct i — Pi) cases (i) and (iii) are identical. 

Equation (5), together with (9)-( 11), relates the shape of the 
time-integrated spectrum to the five parameters of the instanta¬ 
neous spectra and their time evolution, cq, Pi, To, Eo max , and 
f^O.min • 

The general behavior of the result (9)-( 11) is captured well in 
an approximation where the instantaneous spectra are described 
as exponentially cut-off power laws, f y (y) = [A\^ l y ai e~ v , 
where A = T(cq + 1). This corresponds to a Band func¬ 
tion without a high energy power law, i.e., to the case with 
Pi = —oo. The result of the integration using this approxi¬ 
mation has the advantage of resulting in the very simple an¬ 
alytical expression. Equation (10), valid at all energies. The 
effects on the time-integrated spectrum due solely to a finite 
high energy cut-off energy, A’ p i tlnax , is easily seen in (10). By 
letting -Epk.min = 0 and thus y max = oo, the first term on 
the right-hand-side of Equation (10) equals unity. Similarly, a 
non-zero A' p k tllitl will give rise to deviations from the A’ -1 - law. 
Asymptotically, at t/ max <C 1, i.e., at E <C A' p i < ttlitl , the time- 


A/e (E) 


$o ( E \ 

£'0,minr(o;i + 2) \Eq 

,min / 


( 12 ) 


i.e.. the time-integrated spectrum takes the low energy shape of 
the instantaneous spectrum. 

The result of integrating instantaneous spectra described by 
Band functions, i.e.. Equations (9)-( 11), is shown in Figure 1. 
The spectra are in the form of the correction factor, 1(E), and 
are plotted against E/E p k, m ax for different values of rq. The 
time-integrated specific flux spectra, EJVe(E), would be a con¬ 
stant = $ 0 in the delta function approximation, i.e., 1(E) = 1. 
Here, with a more realistic description of the instantaneous 
spectra, there are deviations from the constant To in a bound¬ 
ary interval below A’ p | ! rna x- The width of this boundary interval 
depends on the shape of the instantaneous spectra themselves. 
With an instantaneous photon index of cq = —1, the integrated 
spectrum will never reach the asymptotic level of To. The solid 
curves in Figure 1 show the effect on the To /-£?-spectrum, due 
to the existence of a maximal peak energy. The dotted curves 
also include the effect of a low energy bound to the peak energy. 

As seen in Figure 1, the existence of a maximum and min¬ 
imum peak energy makes the time-integrated photon Band- 
index, ati, that one would obtain when fitting the spectrum, 
take a range of values enclosing the “universal” value of -1. 
This would correspond to part of the “diversity” of ct t i seen by 
Band et al. (1993). 

It is useful to introduce E p (t) as the peak photon energy of 
the instantaneous ENe(E, t) spectra as the existence of such a 
peak corresponds to the assumptions of the presented scenario 
(i.e., cq > —1 and 8\ < —1). In the case of the Band function, 
A/pk and E p are related through E p k = E p (2 + ai)/(l + ap. 
The exponential decay found by LK96 holds for both A’ p and 
E p k as they used a constant cq in their study. 

The total time-integrated specific flux (keV cm -2 ) of a pulse. 
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Fig. 2.— Flux spectrum of the first pulse in GRB 910807 (467) and the best fit model (solid curve); cf. Fig. 1. See the text for the fit¬ 
ting parameters. The data were rebinned to 10 a. The dotted line shows the asymptotic a: t i = —1 power law, with the corresponding 
fitted normalization, <l>o. Note that $o can be directly read off the ordinate. The dashed curve shows the time-integrated spectrum 
for the same parameters as the best fit model but with 7? p k,min = 0. 


T = / Q °° EMe{E) dE, is dominated by the $0 photons cm -2 
in a logarithmic interval near E pru:ix , and we thus expect 
T ~ $o^p,max- Indeed, performing the analytical integrations 
over E, using JVe{E) for the delta function case and the expo¬ 
nentially cut-off power law case, both give T = 4>o(£ ; Pj max — 
E Pt min) ~ ^o^p.max- A similar result is obtained for the Band 
function case when the time-integrated /3 t i <C -2, but as d x \ 
approaches -2, the high energy power law makes an increasing 
and non-negligible contribution. The total energy of a pulse 
coming from an isotropic GRB at a distance d is approximately 
given by 47ni 2< f > o-E , p ,max. making $0 a useful quantity for esti¬ 
mating global energetics. 

3. MULTI-PULSE SPECTRA OF GRBS 

The derived time-integrated spectrum. Equation (5), is for a 
single pulse within the GRB. Most GRBs, however, consist of 
several overlapping or resolved pulses. If the f? p k,min-values of 
all pulses are below the minimum energy of the observed spec¬ 
trum, then all pulses should have a similar low energy slope 
according to our derivation, but their A’pk.max and $0 may 
differ. The sum of a number of time-integrated pulse spec¬ 
tra will lead to a total spectrum with two spectral “breaks”, 
one at the E p k, ma x of the individual pulse with the smallest 
-Epk.max, and the second at the largest -E'pk.max- The time- 
integrated low energy spectrum of k pulses in a burst is the sum 
of the individual pulse spectra, A/E,tot(-E) = X^k-^E.k(-E') = 
E _1 $o,k^k(£-) which would give a ~ — 1 below the 

smallest £pk,max- This slope may, however, not be observ¬ 
able due to the curvature of the individual pulse spectra (see 
Fig. 1) and as the smallest i? p k,max may be close to or below 
the minimum energy of the observed spectrum. Between the 
smallest and the largest T’pk.max. the spectrum will be a super¬ 
position of individual pulse spectra and the spectral shape will 
be determined by the distributions of $0 and T’pk.max among 
the pulses. If, for instance, <l>o is the same for all pulses and 
their -E p k, m ax is evenly distributed in logarithmic energy, the 


spectrum E 2 Ne(E) will be constant, i.e ., the photon index will 
be ~ —2. Above the largest .Epk.max, the spectrum will show 
the high energy behavior of individual pulses. 

The spectral diversity seen when fitting time-integrated spec¬ 
tra of GRBs (Band et al. 1993) is therefore expected to 
mainly be due to the different statistical distributions of <l>o and 
Tpk.max in different GRBs. The spectra of individual pulses, 
on the other hand, are expected to have a spectral shape sim¬ 
ilar to the ones shown in Figure 1. Simple GRBs with few 
pulses may have bumpy spectra (possibly producing previously 
studied line features, see, e.g., Pendleton et al. 1994b), while 
complex GRBs with many pulses will have smooth spectra. 

4. GRB 910807 

The spectra plotted in Figure 1 resemble real GRB pulse 
spectra. To illustrate this and the discussion above, we wish 
to study a bright burst, exhibiting a smooth, long, single pulse, 
as well separated as possible from other pulses. Generally, the 
first pulse in a burst of several separated pulses should be least 
contaminated by soft contributions from previous pulses. There 
are only a few bursts that satisfy these requirements, e.g., the 
first pulse in the bright GRB 910807 (trigger numbers 647; 3rd 
BATSE catalogue, Meegan et al. 1996). Furthermore, this burst 
has very hard time-resolved spectra and no significant evolving 
spectral shape. The spectra only evolve in the sense that the 
peak energy decreases. This motivates the use of our approx¬ 
imation in our deduction of the analytical expressions for the 
time-integrated spectrum, in which we neglected the explicit 
A’pk-dcpcndence. 

In Figure 2, we have plotted the time-integrated spectrum of 
the Large Area Detector (LAD) data for the first pulse of GRB 
910807. It is the specific flux spectrum, i.e., EMe{E), and 
therefore the <f>o value can directly be read off the ordinate. The 
dotted line shows the constant <f>o-spectrum which is then mod¬ 
ified by the /(S)-factor. The solid curve is the best fit of the 
analytical result (5), with /( E ) found by integrating the Band 
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Fig. 3.— The time-evolution of the shape of the instantaneous spectra of the first pulse in GRB 910807 (467), described by the Band 
parameters ct; and A- The dotted line indicates the power law index —1. Note that there is only a weak time dependence. 


function, i.e, equations (9)-( 11) above in §2. The parameters 
of this fit can be compared to the behavior of the time-resolved 
spectra. For instance, their exponential decay, described by the 
decay constant, $o can be compared to the asymptotic normal¬ 
ization of the time-integrated spectrum. 

The best fit to the time-integrated spectrum (see Fig. 2) gives 
the following parameters: $o = 46 ± 18 cm~ 2 , E p k, m ax = 
300 ± 40 keV, E pk min = 110 ± 25 keV, a { = 0.43 ± o’.26 and 
A = -4.2 ± 0.8 with x 2 / d -o.f•= 113.5/113. The an- and $- 
parameters in the fit of the time-integrated spectrum should be 
interpreted as effective values, i.e., the values they would have 
had in the absence of any spectral shape-evolution. 

The time-resolved spectra were fitted with the Band func¬ 
tion with its four parameters: a;, A> E pk , and a normalization. 
The fits were good, having an averaged x 2 /d- 0 .f.= 104/110. 
Figure 3 presents the values of ct; and A of these fits. Note 
that there is no great variation in the shape of the spectra, 
i.e., the a;- and A-parameters remain approximately constant, 
which makes the use of our approximations in §2 relevant. The 
weighted averages of these two parameters are < a x >= 0.42 
and < A >= —2.8. We numerically checked that a weak 
spectral shape-evolution, e.g., from an = 0.6 to 0.2 gives a 
time-integrated spectrum almost identical to that for a constant 
cti = 0.4. In search of T () , we fit the exponential decay law 
(1) to the measured E pk vs. <b, and find <f>o = 50 ± 4 cm -2 . 
E pk varied from 280 ± 20 keV down to 120 ± 12 keV. This 
is merely a factor of two, which results in that the boundary 
effects become important over the whole spectral range. 

Note the expected agreement between the fit of the time- 
integrated spectrum and the parameters from the fits of the time- 
resolved spectra. The study above has shown that by merely 
studying the time-integrated spectrum, information about the 
instantaneous spectra can be extracted. This is made possible 
by knowing in what way the time-integrated spectrum is formed 
from the intrinsic spectra. 


5. EFFECTS OF A TIME-VARYING SPECTRAL SHAPE 

In the discussion above, we assumed that the explicit de¬ 
pendence of fx{x,E pk ) on 7spi,, could be neglected, i.e., the 
shape of the instantaneous spectra is time-independent. This 
allowed us to find illustrative analytical expressions describ¬ 
ing the time-integrated spectrum. The correction to the general 
time-integrated spectral shape, due to such a time dependence, 
is not easily found analytically. It should be pointed out that 
a; has been observed to vary with time (see, e.g., Crider et al. 
1997a). It is, however, obvious that the effects on the time- 
integrated spectra should be small for moderately varying pulse 
shapes, especially in the case of the instantaneous spectra be¬ 
ing hard. A comprehensive study must be made numerically, in 
which the detailed time-evolution of the spectral shape is mod¬ 
eled. Preliminary results show that the analytical expressions 
found above, i.e., the “averaged” spectrum, differ from the nu¬ 
merically integrated spectra, in which a\ is allowed to vary from 
its minimal to its maximal value, by 20 — 30 per cent or less in 
the most cases. To explore how this affects the fitted param¬ 
eters, we have studied two bursts in which ch evolves notably 
and, in particular, approaches -1. It is also of interest to study 
bursts, in which the peak energy varies over a larger range, and 
thus revealing more of the underlying E~ x - behavior, i.e., the 
asymptotic normalization level. 

In the first pulse of GRB 920525 (trigger numbers 1625; 3rd 
BATSE catalogue, Meegan et al. 1996), cti varies between 
~ —0.3 and ~ —1.0, E pk decays from 1250 ± 300 keV to 
150 ± 70 keV, and the To-value found from the time-resolved 
evolution of E pk is 33 ± 1 cm" 2 . For details of the spectral 
evolution, see, e.g., Figure 5 in Ford et al. (1995), Figure 3 
in LK96, and Figure 2 in Crider et al. (1997a). The fact that 
the instantaneous spectra are soft and evolve requires caution to 
be taken when the analytical models presented above are used. 
In minimizing the x 2 of the fit to the data, two equally good 
minima are found (x 2 /d.o.f.= 97/112), both having some dif- 


1 LK96 found ~ 20 cm 2 using LAD5, while we use LAD4, which has a stronger signal. We also use a wider energy range. This illustrates that the <t>o results 
are dependent on the details of the analysis and therefore should be treated with some caution. 
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Fig. 4.— Same as Fig.2., but for the first pulse in GRB 920525 (1625). The large energy interval over which E-p^ evolves allows the 
time-integrated spectrum to reach the asymptotic level of $ 0 



Fig. 5.— Same as Fig.3., but for the first pulse in GRB 921207 (2083). At the end of the studied time-interval, a\ approaches — 1. 



Fig. 6.— Same as Fig.2., but for the first pulse in GRB 921207 (2083) 
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Acuities in describing the actual time-evolution. The first fit 
has <l)o = 27 ± 2 cm" 2 but the intrinsic spectra are found to 
be too hard; ch = +0.6 ± 0.5, and the _E p k-values are too 
small: E v ^ = 800 — 10 keV. The second fit has reasonable 
values for ch = —0.87 ± 0.05 and f3\ = —2.15 ± 0.13, while 
the rest of the parameters, although consistent with the time- 
resolved fits, are poorly constrained: $o = 230 ± 220 cm -2 , 
k'pk.max = 930 ± 330 keV and f? p k,min = 280 ± 140 keV. 
In judging these results, it should be noted that the spectra do 
have significant spectral evolution, which is not accounted for 
in the fitting. Furthermore, the instantaneous spectra are quite 
soft, which makes the constraining of the parameters difficult. 
Figure 4 shows the time-integrated spectrum of the first pulse in 
GRB 920525. The large energy range over which the peak en¬ 
ergy decays and the low value of .E p k,min allows the spectrum 
to reach the asymptotic level of $o- 

Another example, in which a i evolves to soft values is the 
first pulse in GRB 921207 (trigger numbers 2083; 3rd BATSE 
catalogue, Meegan et al. 1996; for details, see Figure 5 in Ford 
et al. 1995, and Figure 1 in Crider et al. 1997c), with a\ varying 
from ~ 0.0 to ~ —0.7 (as shown in Fig. 5), E p k decaying from 
470 ± 25 keV to 105 ± 5 keV and with a <I> 0 -value of 64 ± 2 
cm^ 2 found from fitting the time-resolved exponential decay of 
Apk. As shown in Figure 5 the ct; parameter is constant in the 
beginning of the evolution but approaches ch = — 1 at the end, 
which could cause problems for the constant-a j-approach. The 
time-integrated fit (see Fig. 6), however, correspond quite well 
to the time-resolved results. The fit to the time-integrated spec¬ 
trum gives a; = 0.0+0.2, fa = —3.4±0.2, Spk,max = 505+15 
keV and /:j,k. = 68 + 12 keV, and an asymptotic normaliza¬ 
tion of 4>o = 70 ± 10 cm -2 . The relatively good agreement be¬ 
tween the time-integrated and the time-resolved fits could sug¬ 
gest that the spectral shape actually is constant in time, or at 
least that this is a good approximation, and that the late soft¬ 
ening is caused by erroneous fits due to the low energy power 
law disappearing beyond the observed energy range. Note that 
Crider et al. (1997c) suggest a break in the exponential decay. 
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6. DISCUSSION 

We have presented a scenario in which the exponential de¬ 
cay behavior of the peak energy vs. fluence leads to a time- 
integrated spectrum with an underlying E~ l -behavior. This 
spectral behavior is not necessarily directly visible in a real 
spectrum, as it is modified by various effects, such as the exis¬ 
tence of upper and lower cut-off energies, +j,k.max and f? p k,min- 
We were able to treat these effects analytically. Other effects, 
such as a time-varying spectral shape, require numerical treat¬ 
ment and is left for future work. 

We derived a five parameter expression for the time- 
integrated spectrum. These parameters describe the properties 
of the instantaneous spectra and their time evolution. Fitting a 
time-integrated spectrum with this expression thus gives infor¬ 
mation about the time-resolved spectra. We have illustrated our 
scenario by analyzing spectra from three GRB pulses. In future 
work we will consider a larger sample of pulses. 

Studies of GRB spectra, in particular, those of individual 
pulses are expected eventually to lead to an understanding of the 
physical processes responsible for the GRB emission. We em¬ 
phasize that it is the instantaneous spectra and their time evolu¬ 
tion that reflect these processes, rather than the time-integrated 
spectrum. Unfortunately, most theoretical spectral models as¬ 
sume that it is the time-integrated spectrum that reflects the un¬ 
derlying physical emission mechanism. 
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